Unimpeded tunneling in graphene nanoribbons 
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We studied the Klein paradox in zigzag (ZNR) and anti-zigzag (AZNR) graphene nanoribbons. 
Due to the fact that ZNR (the number of lattice sites across the nanoribbon = N is even) and 
AZNR (N is odd) configurations are indistinguishable when treated by the Dirac equation, we 
supplemented the model with a pseudo-parity operator whose eigenvalues correctly depend on the 
sublattice wavefunctions for the number of carbon atoms across the ribbon, in agreement with the 
tight-binding model. We have shown that the Klein tunneling in zigzag nanoribbons is related 
to conservation of the pseudo-parity rather than pseudo-spin in infinite graphene. The perfect 
transmission in the case of head-on incidence is replaced by perfect transmission at the center of 
the ribbon and the chirality is interpreted as the projection of the pseudo-parity on momentum at 
different corners of the Brillouin zone. 
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Nanoelectronics based on graphene has become a fast growing field [1|, la with a number of technical apphcations 
[3[. Electrons/holes low energy dynamics is described by the Dirac equation valid around the corners of the Brillouin 
zone (Dirac K[K') points) [j-[6|. One of the most striking consequences of the massless nature of the carriers is 
the so-called Klein paradox, which is the unimpeded penetration through high and wide potential barriers pj, 17[ as 
shown in FigiH Graphene provides a convenient medium for experimental verification of the effect. The underlying 
physics of Klein tunneling is based on the notion of chirality as a symmetry between electrons and holes in graphene. 
Formally, each of the branches of the dispersion conical section is characterized by a specific sign of the pseudo-spin 
projection on to the momentum. The Klein effect is related to the conservation of this quantity across the potential 
barrier. 

In this paper, we explore the notion that chirality is necessary for Klein tunneling. For this one must compare 
transport in systems which provide either chiral or non-chiral massless Dirac fermions such as graphene nanoribbons 
(GNR). This are just an infinitely long carbon nanotubes unzipped along either armchair (ANR) or zigzag (ZNR) 
direction and then flattened. The ANR boundary conditions mix the valleys at K {K') points giving the dispersion 
strongly depended on the number of carbon atoms (TV) across the ribbon (modulo 3) [8| producing semi-metallic 
or insulating subbands. The fermions in the ANR are chiral and the potential barrier is perfectly transparent to 
electron scattering with angles close to normal incidence, in agreement with the Klein effect in the graphene sheets. 
This makes ANR a basis for graphene switches, random access memory and GNR field-effect transistors. When the 
electron motion is confined by zigzag edge, surface states describe the lowest energy bands of the ZNR providing the 
right (left) transmission channels [9[[l0| at K {K') as shown in Fig{TJ Transmission through the barrier is possible 
if the gate potential acts as an inter- valley scattering source. Unlike the ANR, the ZNR lowest energy dispersion 
does not show substantial dependence on N. However, the transmission drastically depends on this parameter. 
In the single channel regime, it has been predicted that the transmission is allowed (forbidden) for odd (even) N . 
This calls for referring to the case with odd number of atoms as anti-zigzag nanoribbon (AZNR) This current 
blocking was originally interpreted in analogy with the spin- valve effect in ferromagnetic junctions, thereby gaining 
the name of valley- valve effect [l2|, Il3|. Later, however, it was attributed to symmetric (asymmetric) coupHng of 
states localized at opposite edges (sublattices) of the GNR [l^. Owing to the fact that the model for the ZNR 
and AZNR configurations are indistinguishable at the level of the Dirac equation [l5[ , the latter must be augmented 
with a pseudo-parity operator whose eigenvalues provide the correct dependence (in agreement with the tight-binding 
model) of the sublattice wavefunctions on TV. Bellow, we show that when it comes to Klein tunneling in AZNR, 
pseudo-parity plays the role of the pseudo-spin in infinite graphene. Perfect transmission in case of head-on incidence 
is replaced by perfect transmission at the center of the ribbon. The chirality in AZNR is related to the projection of 
the pseudo-parity on the pseudo-spin around and K'^. 

The energy band structure of GNR depends on the shape of their edges. We shall focus on the ZNR configuration 
since the fermions in this geometry are not chiral [l5|. At {K') points, there is an excess right (left) moving 
conducting channel for 7r(7r'^) bands. In the spirit of single channel regime [10|, these bands, and only them, will be 
considered further. Conventionally, these are given by k • p approximation for the tight-binding model [l| yielding the 

following electron + (hole — ) dispersion valleys: V{x) — E = — k^. Here the energy is measured in units of the 

hopping parameter t = 2.7a/3/2 eV between nearest-neighbor carbon atoms and the wave vectors are in units of the 
inverse distance between them ao = 1.42 A. The hard wall boundary conditions mix the transverse and the longitudinal 
electron momenta kx = ky coth (Wky)^ where the GNR width isW = 2N — 1 for even number N of carbon atoms across 
the ribbon and W = |(A^ — 1), otherwise. The chemical potential provided by the underlying metal contact (insulating 
strip, carbon nanotube) [H, ll^ assumes the form of a sharp n — p — n junction: V{x) = Vq [0-^{x) — Oj^{x — D)]. In 
this notation, Vq is the height of the potential barrier and 9j^{x) is the Heaviside step function. In this way, we obtain 
the following equations: 

As shown in Fig. [TJ this alternating doping provides the scattering {p— doped) region sandwiched in between left 
(right) leads (n— doped). The electric current through the potential barrier is given by the incident electron of energy 
E in the range [0, A], where 2 A is the energy separation between the top of the next highest valence and the bottom 
of the next lowest conduction band (parabolically shaped V, (V"^) bands in Fig. [1]). To stay in the single channel 
regime, where the scattering occurs within lowest conduction tt and highest valence tt"*" bands, we require the barrier 
height Vq in the range \E — Vq\ < /S.. 

The scattering problem with the sharp edge junction is fully determined by the wavefunction for the tt and tt^ bands. 
The corresponding eigenfunctions for the two sublattices (A,B) can be factorized into the element-wise product of 
longitudinal, transverse, pseudo-parity and electron-hole parity two component wavefunctions as: 

^±k,n{x, y)= {^^] = 0±k(^) o x{y) o V±K o s (1) 
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Here, the longitudinal wavefunction depends on the parity Xn of the number of carbon atoms across the ribbon: 

, / A — exp (±i(iir + \ 

cP^kM^) _ -= ( ^^^^^ ^ ^^^^^ ^ ^^^^ 1 (2j 



, where Xeven = and Xodd = a/3/2. Its presence stems from the suitable shape of the unit cells as depicted in 
FiglH The choice of ± defines the right (transmitted) and the left (reflected) movers (direction of the electron group 
velocity for E — V{x) > 0), respectively. The shape of the unit cells and, consequently, the form of the longitudinal 
wavefunction ascertains translational symmetry of the wavefunction ([T]) along the ribbon. The transverse wavefunction 
component is 

y(y) - ( \ ( Si^h {ky{Wl2 + y)\ \ 

^^y' ~ \smh{2Wky) - 2Wky J \ sinh [ky{W/2 - y)] ) ^''> 

,which satisfles the Dirac equation at = ±27r/3A/3 as well as the hard- wall boundary conditions across the ribbon. 
In contrast to inflnite graphene or the ANR, these bands describe localized bound edge states for Wky = and 
delocalized bulk states for 5Re Wky = 0, Wky < tt; while V(V'^) bands correspond to Wky = 0, Wky > tt. 
It is interesting to notice that, in contrast to inflnite graphene, the head-on incidence (^m ky <C kx occurring at 
Wky nir) is prohibited due to mixing of the longitudinal and transverse wave vectors. 
There are two types of electron (hole) states in the ribbons given by the bonding and anti-bonding mixing of the 
sublattice wavefunctions. These strongly depend on the mirror symmetry (or lack of it) of the GNR, which follows 
from the tight-bing model To incorporate this feature into our model, we introduce the concept of pseudo-parity 
similar in nature to that in Ref. [l2j : 

"^^^ " ( {±K3V3/27r)^exp[i^m {Wky)] ) 

The pseudo-parity component of the wavefunction ([T]) determines wether the electron (hole) wavefunctions are either 
even or odd, forming bonding or anti-bonding states. Formally these can be deflned by the following transformation 
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The pseudo-parity ensures alternation of the wavefunction sign for different bands as shown in Fig. [TJ In fact, the 
parity of N determines whether or not the bonding and anti-bonding wavefunctions switch between conduction and 
valence bands at K {K') points. The electron- hole parity is comprised of the eigenvalues of the time-reversal operator: 

sgn(E-F(x)) ) ' 

where positive sign correspond to electrons in the lower conduction tt band, while negative sign stands for the holes 
in the highest valence ir'^ band. 

The problem of tunneling through the n—p — n barrier can be described as follows. An incident particle belonging 
to the right moving channel {K point) for x < strikes the potential barrier at x = 0. The nature of the scattering 
determines the transmission probability density by employing continuity of the wavefunctions ([1]) at the junctions 
and conservation of the pseudo-parity iri±. The reflected wave contributes to the left-moving channel (at K'). It is 
clear that the necessary condition for non-zero transmission is the change of the wavefunction parity under reflection. 
This fails for N even (ZNR conflguration in Fig(T]) and the transmission vanishes identically t = 0. The sublattice 
wavefunctions ([1]) for tt band in the n-doped region and tt"*" band in the p-doped region have opposite sign for the 
pseudo-parity (opposite bonding) [14]. The situation is drastically changed for AZNR (A^-odd) conflguration. The 
reflected wavefunction changes its pseudo-parity while the transmitted one keeps the parity of the incident particle 
(tt band at K and tt^ at K' have the same pseudo-parity). The transmission probability density is given by 

t{y) = 7 — ^ T — (7) 

Here, we introduced P = exjp [iitDXn / \fS) and A = sb^iSb,2Xa,iXa,2Xb,iXb,2] the sub-indices A^B on the 
wavefunction components denote the corresponding sublattices. The gated {p doped) and n doped regions are specifled 
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by the sub- indices 2 and 1 respectively, for instance kx,2^ stand for the wave vectors along the ribbon in the p and 
n doped regions. The transmission probability (and hence, in the single-channel regime, the conductance in units 

w 

of 2e^/h) is given by: \t{kx^i^kx^2)\^ ~ W S ^V'^iy)'^^ iv) shown in Figl2l It jumps up to unity at 5Re ky^2D = nn 



for any integer n in case of surface bounded holes, and at ky^2D = n7r/2 for the bulk-like holes. For the latter, 
there are additional resonances at ky^2D ^ nn. This situation corresponds to the quantum-dot-like states in the 
gated region [l8|. These resonances are similar to those in infinite graphene as shown in Figl3l The transmission is 
also ideal on the diagonal kx^i = kx,2 and corresponds to the absence of the gating potential. 

Above we analytically demonstrated that when the Dirac fermion picture is supplemented with the concept of 
the pseudo-parity for ZNR (AZNR) it removes the ambiguities of the ballistic transport and correctly describes the 
valley- valve effect (current blocking). We now turn to the role of pseudo-parity in the Klein effect in nanoribbons. As 
mentioned before, the Klein paradox in infinite graphene or metallic ANR is related to the conservation of pseudo- 
spin projection on the particle momentum, quantity known as chirality. That is the through the barrier unimpeded 
transmission occurs at the head-on incidence as indicated in Figl3](a) for kx^iD = 27r for all values of the gate potential 
Vb (any kx^2D). Transformation from the sublattice basis to the bonding/anti-bonding representation yields that the 
head-on incident electron is in the purely bonding state (^^ = 0) However, the fermions in tt and n"^ bands for 
ZNR and AZNR do not satisfy the conventional (pseudo-spin based) definition of chirality, thus often refereed to as 
non-chiral. Furthermore, the net transmission does not demonstrate such an effect. We now turn our attention to the 
transmission density given by Eq. ([7]). It shows that at the center of the nanoribbon, transmission density is always 
perfect t{y = W/2) = 1, regardless of the gate potential Vq. This effect is intimately related to the Klein paradox in 
infinite graphene. Since, at the nanoribbon center, the incoming wavefunctions from the sublattices A and B form 
the bonding electron state (^a(^ = W/2) = 0). This holds true for both edge-bound and bulk states of the tt and 
7r~^ bands. It follows that we can relate the unimpeded transition density to conservation of the pseudo-parity rather 
than conservation of the pseudo-spin in infinite graphene. In AZNR, the pseudo-parity itself defines the chirality of 
the fermions. The valley- valve effect of the current blockage is the other manifestation of the Klein paradox through 
the pseudo-parity defined chirality in ZNR. 

It was also shown that with the use of density functional theor y, it is possible to obtaina more accurate band 
structure, which shows a gap in the ground state of zigzag nanoribon [l^. It was also shown by Son et al that the 
gap is inversely proportional to the width of the ribbon. So, for a ribbon which is wide enough the bandgap becomes 
negligibly small. In our model calculations, the pseudo-parity of the nanoribbon is more significant compared with 
intrinsic spin. To keep the model simple and to focus on the leading effect we have neglected spin altogether. 

In conclusion, the Klein paradox is unimpeded tunneling of the purely bonded Dirac electron state across arbitrary 
wide gated region in infinite graphene. Its another manifestation is perfect refiection in the graphene stacks. The 
paradox is conventionally interpreted in terms of the chirality conservation. The latter is defined as a projection of 
pseudo-spin on the direction of motion s • p for the two branches originated from the sublattices. In more general 
sense, the term chirality is often used to refer to an additional built-in symmetry between electron and hole parts of 
the spectrum. However, in chemistry, the term chiral is used to describe an object that is non-superposable on its 
mirror image. These two definitions cannot be related to infinite graphene since this system has no mirror refiection 
symmetry. The situation is not the same for the GNR. In accordance with the chemistry definition, ZNR's are 
non-chiral (the mirror plane is perpendicular to the GNR plane). On the other hand, AZNR is a chiral object in 
the chemistry sense. It is also chiral in the sense of pseudo-parity (two opposite projections of the pseudo-parity 77 
on the pseudo-spin s). As we have demonstrated, its conservation across the n — p — n junction defines unimpeded 
electron tunneling at the AZNR center. Its conservation under the refiection is a source of valley- valve effect (perfect 
refiection) in ZNR. Both are facets of Klein paradox in nanoribbons. 
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I. APPENDIX 



In this Appendix, we consider the same problem for armchair nanoribbons. The wave function has the following 
form: 

^ = ( "^^^ = ( ■ ^) exp(z/c^x)0A(r) + exp(zK' • r) exp(i/c^x)^^(r) \ 

J \ exp(zK • r) exp{ikxx)(j)B{r) + exp(zK' • r) exp(z/c^x)^'^(r) J ^ ^ 
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Here 

( \ exp{-i(l)) exp{ikny) \ f ^aA^^/^ - exp{-i(l)) exp{-ikny) \ 

J \ exp{ikny) J J \ exp{-ikny) J ^ ^ 

The wave function ^ must be continuous at x = and x = D corresponding to the boundaries of the region with 
potential V{x). 

^(0, y; kx, ky] si; 0) + r^(0, y] -k^, ky] 5i; tt - 0) = a^(0, y] Qx, ky] S2] 6>) + 

+6^(0, y] -Qx, ky] si; TT - 6) 
a^{D, y; Qx, ky; 52; 0) + b'^{D, y; -Qx, ky] 5i; tt - 0) = t^{D, y; kx, ky; 5i; 0) (10) 

In the armcair case we have the same expressions connecting kx^ kn-, and (j): 

kx = kfCOs{(j)) ; kn = kfsin{(j)) 

= ^{E - VoY - kl ; ^ = tan{d) (11) 

Also the phase for reflected waves is tt — (/) compared to the phase (j) for incoming wave beacuse only this change of 
phase will ensure opposite kx with the same kn accourding to equataions (pTj) 

In the simplest approximation when there is no K and K' valley mixing, the above wave-function reduces to the 
four component vector 

/ exp(-z^)exp(i/cn^) 
T / \ s exr>(—iknV) 

^ ' I -exp(-z0)exp(-z/Cn^) 
V SQxp{-ikny) 

This gives rise to eight equations instead of the usual four which we need for obtaing the unknowns r, a, 6, t. Our 
equations may be divided into four pairs, for which both equations in each pair are equivalent. Consequently, we 
have the same set of simultaneous equations to solve as for infinite graphene case, except kn accross the nanoribbon 
is quantized and given by 




- 3ao • ^^^^ 

A straightforward calculation leads to the coefficient for transmission i = t{(t)^ 6) is and hence the transmission 
probability given by pLi] 

T = tH= co.^(<^) cos\e) _ ^^^^ 

[cos{Dqx)cos{(j))cos(0)] + siirfi{Dqx) [1 — siS2sin{(j))sin{0)\ 

In the limit of a high potential barrier, i.e., \Vq\ ^ the transmission (fT4|) will becomes 



1 — cos^{Dqx)siii?{(j)) 

The case of special interest is = (3A/' + l)ao which corresponds to the metallic subbands structure - there is no gap 
between the highest valence and the lowest conduction subbands. In this case /cn = and as a result = and the 
expression (p!5]) will give T = 1 - perfect transmission. For non-metallic case we have sin{(j)) = 7r/kfW. 
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(c) 



FIG. 1. (Color online). Panel (a) shows infinite graphene n — p — n junction and corresponding chiral dispersion curves. The 
TT, tt'^ bands are depicted by black lines; the nearby parabolic shaped (blue) curves denote V, V"*" bands. The different values 
of the pseudo-spin projected on the momentum (eigenvalues of cr • p given by sbPx/\p\) are shown in solid and dashed 
curves. Panel (b) describes 9-AZNR n — p — n junction and chiral particle dispersion. The two different values of pseudo-parity 
projection on the pseudo-spin 1 — s"^ • 77 are presented in solid (dashed) curves. Panel (c) shows 8-ZNR configuration and 
non-chiral dispersion curves. For all panels dashed rectangles represent the unit cells. 
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FIG. 2. (Color online). The transmission |t(/cy,i, /cy,2)|^ for 9-AZNR non-chiral tt electron. Panels (a.l)-(d.l) correspond to 
over the barrier transmission E > Vq. Panels (a.2)-(d.2) represent through the barrier E < Vo transmission, A^y,i(2) implies 
that the electron (hole) is in the bulk state, while A^y,i(2) signifies the surface-bound states. 
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FIG. 3. (Color online). The transmission \t{kx,i, kx,2)\'^ for an infinite graphene. Recall that conventionally the transmission 
is given in terms of the incident and refraction angles (01, ^2) between the transverse and longitudinal wave vectors in n- and 
p-doped regions correspondingly. The conversion is simply given by Xkx^i = 27rsin((/)i) cot 02, Xkx,2 — 27rcos(/)i. As in Ref.[i|], 
the incoming electron wavelength is taken to be equal to the width of the p-doped region. Through the barrier E <Vo (a) and 
over the barrier E >Vq (b) transmission of the chiral electron. 



